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Aims. The cluster 1ES0657-556 is an ideal astrophysical laboratory to study the distribution and the nature of Dark Matter 
because this last component is spatially separated from the intrac luster gas. We show that microwave observations can provide 
crucial probes of Dark Matter in this system. 

Methods. We calculate the expected SZ effect from Dark Matter annihilation in the main mass concentrations of the cluster 
1ES0657-556, and we estimate the sources of contamination, confusion and bias to asses its significance. 

Results. We find that SZ observations at v w 223 GHz can resolve both spatially and spectrally the SZdm signal and isolate it 
from the other SZ signals, and mainly from the thermal SZ effect which is null at v ~ 220 — 223 GHz for the case of 1ES0657-556 . 
We conclude that SZ observations with J$ arcmin resolution and ^ fiK sensitivity of 1ES0657-556 are crucial, and maybe 
unique, to find direct astrophysical probes of the existence and of the nature of Dark Matter, or to set strong experimental limits. 



1. Introduction 

Dark Matter (DM) annihilations in the halo of galaxies 
and galaxy clusters have crucial astrophysical implica- 
tions. In fact, if DM is constituted by weakly interact- 
ing massive particles (for which the leading candidate is 
the lightest supersymmetric particle, plausibly the neu- 
tralino x) •> their annihilation produces secondary particles 
(e.g., neutral and charged pions, secondary electrons and 
protons, neutrinos) that give rise to various astrophysi- 
cal signals. These are, among others, observable fluxes of 
positrons, antiprotons, gamma rays, neutrinos, as well as 
signals due to secondary electrons which cover the whole 
e.m. spectrum (see, e.g., Colafrancesco et al. 2006 for 
details): synchrotron radio emission (in the intra-cluster 
magnetic field), bremsstrahlung emission (if there is co- 
spatial intra-cluster gas), inverse Compton emission due to 
the up-scattering of CMB photons and hence a specific SZ 
effect (as first noticed and derived by Colafrancesco 2004). 
The spatial and spectral intensity of the astrophysical sig- 
nals coming from XX annihilation is expected, however, to 
be confused or even overcome by other astrophysical sig- 
nals originating from the intracluster (IC) gas and/or from 
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the relativistic plasmas present in the cluster atmospheres, 
especially when all these components are co-spatially dis- 
tributed with the DM component. This situation occurs 
in most galaxy clusters (see discussion by Colafrancesco 
et al. 2006 for the case of Coma). 

An ideal system to detect DM annihilation signals would, 
therefore, be a cluster with a clear spatial separation be- 
tween the various matter components. This is, indeed, the 
case of the cluster 1ES0657-556 where the spatial dis- 
tribution of DM is clearly offset w.r.t. that of the IC 
gas (Clowe et al. 2006). The two baryonic clumps of 
hot gas emit X-rays by (thermal) bremsstrahlung, as ob- 
served by Chandra (Markevitch et al. 2002, 2004). The 
shock observed in the western-most region of the clus- 
ter (Markevitch et al. 2002) might be the site of high 
energy emission from particles accelerated at the shock. 
Hard X-ray emission from the direction of 1ES0657-556 
has been marginally detected by Rossi-XTE (Petrosian et 
al. 2006) but its angular resolution is not sufficient to give 
any information on the spatial distribution of this emis- 
sion. No gamma-ray emission has been detected from this 
system with EGRET. The extended radio halo associated 
to this cluster (Liang et al. 2000) has a surface brightness 
slightly elongated along the direction of the two X-ray 
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clumps, but without clear signatures of radio-brightness 
enhancements at the DM clump locations. Finally, the SZ 
maps of 1ES0657-556 obtained with ACBAR (with - 4.5 
arcmin FWHM resolution, Gomez et al. 2003) are quite 
smooth and regular with no evidence of enhancement at 
both X-ray and/or DM clump locations. In this Letter, 
we will compute the specific feature of the SZ effect (here- 
fater SZE) produced by DM annihilation, SZdm, in the 
cluster 1ES0657-556 and we will show that it is possi- 
ble to detect such SZdm signal with a specific observa- 
tional strategy. The relevant physical quantities are cal- 
culated using H = 70 km s" 1 Mpc -1 and a flat, ACDM 
(O m = 0.3, = 0.7) cosmological model. 

2. The complex SZ effect in 1ES0657-556 

The various SZ signals expected from the subsystems of 
the cluster 1ES0657-556 are: i) the SZdm effect, which 
is expected to be located at the two DM clumps; ii) the 
thermal SZ effect (SZ t h) which is expected to be located 
at the two X-ray clumps. We will compute in the fol- 
lowing these two sources of SZE and we will also discuss 
the possible sources of contamination, bias and confusion. 
The general expression for the SZE which is valid in the 
Thomson limit for a generic electron population in the 
relativistic limit and includes also the effects of multiple 
scatterings and the combination with other electron pop- 
ulation in the cluster atmospheres has been derived by 
Colafrancesco et al. (2003). This approach is the one that 
will be used for the derivation of the SZdm effect induced 
by the secondary electrons produced by XX annihilation 
(see the original derivation by Colafrancesco 2004) as well 
as for the SZ^ produced by the hot X-ray emitting gas. 
According to these results, the spectral distortion of the 
CMB spectrum induced by a population of electron with 
momentum distribution f e (p) can be written as 



Ai( X) = 2 ^g£ a m , 

where To is the CMB temperature and 
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in terms of the pressure P e contributed by the specific 
electron population. The spectral function g(x), with x = 
hv/k&To, can be written as 



m e c" 



io(xe s )P(s)ds — io(x) 
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in terms of the photon redistribution function P(s) and 
of i (x) = 2(k B T ) 3 /(hc) 2 • x 3 /(e x - 1). Here r = J d£n e 
is the optical depth of the electrons with number den- 
sity n e , and (e) = ^ J P e d£ = / °° dpf e {p)\pv{p)m e c is 
the average energy of the electronic plasma. The pho- 
ton redistribution function P(s) = J dpf e (p)P s (s;p) with 
s = ln(i//z/), in terms of the CMB photon frequency in- 
crease factor v' jv, depends on the electron momentum 



distribution / e (p), where the momentum p is normalized 
to m e c. The CMB temperature change produced by the 
SZE is finally given by 



(4) 
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The specific SZdm and SZ t h effects for the various elec- 
tronic components in the cluster 1ES0657-556 are com- 
puted following the approach previously described. 

2.1. The SZ DM effect in the cluster 1ES0657-556 

The calculation of the secondary electron spectrum from 
XX annihilation in galaxy clusters has been already pre- 
sented in details by Colafrancesco & Mele (2001) and 
Colafrancesco et al. (2006), and here we will only recall 
the relevant steps necessary for the present purpouses. 
We assume, for simplicity, a spherical DM halo model for 
each DM clump of the cluster 1ES0657-556 , as indicated 
by the lensing maps derived by Clowe et al. (2006), with 
DM density profile given by g(x) — x _r? (l + x) 7 ? - ^, with 
x = r/r s . Values r\ = 1 and £ = 3 reproduce the Navarro, 
Frenk & White (1997) density profile. 
The neutralino number density profiles n x (E,r) = 
n x $(E)g{r) of the two DM clumps have been calculated 
following the approach described in Colafrancesco et al. 
(2006) with a NFW DM density profile and the following 
structure parameters: M vir = 1O 15 M , R v i r = 1.97 Mpc 
and c v i r — R V i r /r s =5.66 (for the larger East DM clump); 
M vir = 6.25 • 1O 13 M , R vir = 0.784 Mpc and c vir = 7.56 
(for the "bullet" West DM clump). 

The x annihilation rate in the DM clumps is R = 
n x (r)(av)o , where (ctv)q is the XX annihilation cross 
section averaged over a thermal velocity distribution at 
freeze-out temperature. The range of neutralino masses 
and pair annihilation cross sections in the most general 
supersymmetric DM setup is extremely wide (see dis- 
cussion in Colafrancesco et al. 2006, 2007). We consider 
here, specifically, the neutralino models worked out in 
Colafrancesco et al. (2006) with M x = 20,40 and 81 GeV 
and with their specific values of (<jv)o- 
The electron source functions Q e (E,r) oc (av) qti x {E , r) 
for the specific neutralino model considered here have been 
derived in Colafrancesco et al. (2006) and the time evolu- 
tion of the electron spectrum is given by the equation 



dn e (E,r) 
dt 



_d_ 

dE 



n e (E,r)b(E) 



Qe(E,r), 



(5) 



where spatial diffusion can be safely neglected in cluster- 
size DM clumps (Colafrancesco et al. 2006). The function 



+b°couinth (1 + Iog(7/nt h )/75) 
+bb rem nth (log( 7 /n th ) + 0.36) 



(6) 



gives the energy loss per unit time at energy E where n t h 
is the mean number density of thermal electrons in cm -3 , 
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7 = E/m e c 2 and b° IC ~ 0.25, b° ayn ~ 0.0254, b° Coul ~ 6.13, 

6 brem - L51 > a11 in units of 10 ~ 16 GeVs" 1 . The equi- 
librium spectrum n e (E,r) obtained solving eq.(0) allows 
to calculate the SZdm effect. FigJT] shows the CMB tem- 
perature change AT evaluated at the centers of the two 
DM clumps for different values of M x . The SZ dm signals 
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Fig. 1. The CMB temperature decrement produced by the 
East (left panels) and West (right panels) clumps of gas (red 
shaded areas which take into account the effect of IC gas tem- 
perature uncertainties) and DM with neutralino mass M x = 20 
GeV (green), M x = 40 GeV (blue) and M x = 81 GeV (cyan). 
The lower panels show a zoom of the SZdm effect only. 

are overcome by the SZ^ signals at low (y ^ 200 GHz) 
and high (y ^ 230 GHz) frequencies, while they dominate 
in the frequency range where the zero of SZ^ is found, 
i.e. at v ~ 223 GHz for the East clump and at v « 219 
GHz for the West clump. At such frequencies the SZdm 
temperature decrement takes values « —21.1,-10.6,-0.3 
fjK (East clump) and values « -8.4, -4.2, -0.1 /iK (West 
clump) for M x = 20,40,81 GeV, respectively. 

2.2. The SZ th effect in the cluster 1ES0657-556 

The SZth from the two X-ray clumps of the cluster 
1ES0657-556 is also computed from the general approach 



delineated in eqs.(l-4). Following Markevitch et al. (2002, 
2004) and Tucker et al. (1998), we use the simplifying as- 
sumption that the IC gas distribution can be described 
by isothermal spheres fitted by a /^-profile. For the East 
X-ray clump we adopt kT e = 14 keV (Markevitch et al. 
2002, 2004) and a ^-profile with r c = 257 kpc, f3 = 0.62 
and a central IC gas density n t h,o = 6.71 • 10 -3 cm -3 
(Tucker et al. 1998). For the West X-ray clump we adopt 
kT e = 6 keV and a /3-profile with r c = 125 kpc, /3 = 0.7 
and a central IC gas density n t h,o = 2 • 10 -2 cm -3 , ne- 
glecting the additional central power-law density profile 
induced by the presence of the shock (Markevitch et al. 
2002, 2004). Thus, SZ th from this last X-ray clump should 
be considered as slightly underestimated. 
The possible non-equilibrium and non-thermal effects in- 
duced by the merging shock are not crucial for the aims of 
this paper since the shock is located outside the two X-ray 
clumps and at a large angular distance from the two DM 
clumps (see Markevitch et al. 2002, 2004). Possible non- 
thermal SZE produced by shocks is, however, expected to 
be at most at a level of a few /iK (see Colafrancesco et al. 
2003), and therefore it could be, at most, a source of bias 
for the thermal SZE at the location of the X-ray clumps. 
The SZ t h temperature change for the two X-ray clumps is 
also shown in FigHJ A major difference between the SZ^ 
and SZdm spectral functions is the different position of 
the zero of the SZE which is moved to higher frequencies 
in the case of the DM produced electrons with respect to 
the case of the thermal distribution. As a consequence, 
the SZdm effect appears as a negative contribution to the 
overall SZE at all the frequencies which are relevant for 
the SZ experiments, x ~ 0.5 — 10. 

2.3. Contamination, bias and confusion 

Possible sources of confusion for the microwave search 
for DM in 1ES0657-556 are: CMB anisotropy, emission 
of the IC medium, emission of galaxies and AGNs along 
the line of sight, thermal and kinematic SZE from the 
X-ray clumps. For the estimates of the foreground con- 
tamination we consider as peaks of the DM distribution 
producing SZE the regions included in the n = 0.3 contour 
of Fig.l in Clowe et al. (2006). These are roughly 1.3' in 
diameter, and are separated by ~ 4.8'. 
The measurement we propose must be carried out with 
an instrument with angular resolution of ~ V FWHM, 
able to resolve the two SZdm peaks, and to distinguish 
them from the thermal and kinematic SZE of the two X- 
ray clumps. Larger scales are not significant for the mea- 
surement, and must be removed either in the measure- 
ment or in the data analysis. This will remove most of 
the CMB anisotropy. In fact, integrating over the power 
spectrum of CMB anisotropy best fitting current measure- 
ments (including lensing effects), we find that an angular 
high-pass filter removing scales ^ 3' will cut multipoles 
£ ^ 3600, leaving a residual rms fluctuation of primary 
CMB anisotropy of 0.8 /iK rms. This is further reduced 
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by the finite angular resolution of the telescope. 
In the two SZdm peaks we find 14 galaxies in the area of 
the East peak and 6 galaxies in the area of the West peak 
(see Barrena et al. 2002). To obtain an estimate of the 
mm- wave emission of these galaxies, we take the mm- wave 
spectrum of an "average" , "normal" galaxy (like M99) and 
redshift it to z=0.3 . Then we compute the flux in the 
bands centered at 90, 145, 217, 270, 345, 545 GHz and 
we find the total flux by multiplying it by the number of 
galaxies counted in the West and East SZdm peaks. In the 
West peak, we obtain a total flux of 0.2, 0.4, 1.0, 1.6, 3.4, 
19 mJy, respectively, in the bands quoted above. These 
signals correspond to a CMB anisotropy of 4.0, 4.8, 8.0, 
13, 40, and 760 \iK for a 1' FWHM beam. These are com- 
parable in size to expected SZdm signal, at least in the 90, 
145, and 217 GHz bands, but are positive signals, while 
the SZdm signal is negative at the same frequencies. The 
same holds for AGNs: there is only one known radio source 
(SUMSS J065837-555718) listed in the NED database in 
the area of the East peak, and none in the area of the West 
peak. This source has been observed with the MOST ra- 
dio telescope, and features a power law SED, with slope 
a ~ —0.9. If the same SED is used to extrapolate to the 
frequencies of interest here, the total flux is safely negli- 
gible with respect to the expected SZdm in all the bands 
considered here: for a 1' FWHM beam the signal is below 
1 [iK in the 145, 217, 270 and 345 GHz bands. 
Unresolved background galaxies produce a flux with posi- 
tive and negative fluctuations with respect to the average. 
For a 1' beam at 220 GHz, we get a rms signal of ~ 
estimated using the model of Lagache et al. (2004). 
The SZ t h emission from the IC gas is offset with respect 
to the DM peaks, so that only the tail of the signal is 
present in the SZdm peaks. The SZ^ emission is cer- 
tainly affected by the complex temperature distribution 
of the X-ray clumps (see Markevitch et al. 2002) produc- 
ing fluctuations in the SZ maps (since AT ex j d£n e T e ). 
This makes the maps at low (150 GHz) and high (350 
GHz) frequency shown in Fig. 2 as only indicative of the 
observable SZ signal. However, we stress that the SZ maps 
at v « 219 — 223 GHz are quite realistic because at 223 
(219) GHz for the East (West) X-ray peak, the SZ th is 
null, and temperature uncertainties do not sensitively al- 
ter this result (see FigJT]): this is, henceforth, the privileged 
band for this measurement. 

The (AT /To) kin at the clump centers is Si 2.7 • 
10- 6 (r/10- 3 )(l/ p /10 3 A:m < s- 1 ) and is maximum at 217 
GHz. Small values of AT^ n are found at the location of 
the two DM clumps because: i) the optical depth of the 
secondary electrons is small ^5 a few 10 -5 ; ii) the optical 
depth of the residual IC gas at the DM peak locations is 
small a few 10 -4 , and hi) the peculiar velocity of the 
X-ray clumps and DM clumps along the line of sight is rel- 
atively small (Markevitch 2004), since the major merging 
event occurs almost completely on the plane of the sky, i.e. 
mostly perpendicular to the line of sight. Even an extreme 
value V p « 600 km/s (well above the rms value « 300 ±80 
km/s found for galaxy clusters, see Giovanelli et al. 1998) 



would produce a total AX/^ n ~ 0.4 /iK (~ 0.2 /iK) in the 
East (West) DM clump. In conclusion, being the two DM 
peaks offset with respect to the SZ^ peaks, the kinetic SZ 
effect is at most a marginal residual signal at the location 
of the two DM peaks of the cluster 1ES0657. 
In these conditions the detectability of SZdm is limited 
mainly by the fluctuations of the unresolved galaxies back- 
ground. Negative or positive fluctuations can be found 
at the locations of the SZdm peaks, affecting its detec- 
tion. From realistic scan simulations at 220 GHz with a 1' 
FWHM beam, we estimate that a 2-a (3-cr) detection has 
a likelihood ^ 95% if the SZdm in the beam is 6.5/iK 
( Si 10/iK). These limits are comparable to the expected 
SZdm signal estimated in Sect 2.1. Complementary obser- 
vations at higher and lower frequencies should confirm the 
nature of the signal, which in the case of the unresolved 
galaxies increases strongly with frequency, while is basi- 
cally constant for the SZdm There are several instruments 
expected to provide SZ surveys from the southern hemi- 
sphere in the near future. The 10-m South Pole Telescope 
(SPT, Ruhl et al. 2004) will exploit the excellent envi- 
ronmental conditions of the antarctic winter to produce 
deep surveys of SZ clusters at 90, 150, 220 and 270 GHz. 
With a FWHM of ^ 1' this instrument is ideally suited 
to perform the DM search outlined here. The sensitivity 
of the SZ survey will be ~ 10/aK for each V pixel at 150 
GHz: we suggest deeper observations for 1E0657-558. A 
sensitivity of ~ 0.5/iif per V pixel at 150 GHz can be ob- 
tained in 4 days of integration. Increasing atmospheric and 
bolometer noise prevents reaching the same sensitivity in 
the 220 and 270 GHz bands (see also Gomez et al. 2003, 
Runyan et al. 2003 for the case of experiments working at 
larger angular scales). Still, the sensitivity should be good 
enough to see the two ~ 10/nK negative signals of SZdm- 
Similar considerations apply to the Atacama Cosmology 
Telescope (Kosowsky 2003), to the APEX bolometer ar- 
rays (Lee et al. 2004) and to the BOLOCAM instrument at 
the CSO (Glenn et al. 1998). Fig |2] shows the CMB maps 
expected from an observation of 1ES0657-556 at v = 150, 
223 and 350 GHz (with bandwidths of w 25%) with the 
SPT. The SZdm effect clearly emerges localized at the 
two DM clumps for v = 223 GHz where SZ^ vanishes 
for the brighter X-ray clump with kT e — 14 keV. This is 
the optimal frequency to detect the SZdm from the clus- 
ter 1ES0657-556 . Such an experiment is able to set strong 
constraints on DM particle mass in the range M x ~ 20 — 50 
GeV, while the SZ signals for M x ^ 80 GeV is unobserv- 
able. Higher frequency measurements are needed to mon- 
itor foreground contamination. Observations of 1E0657- 
558 with the OLIMPO experiment (Masi et al. 2005), 
using bolometer arrays at 145, 217, 345 and 545 GHz, 
with beam FWHM of 3.3', 2.2', 2.2', 2.2' respectively, will 
be able to resolve the two SZE peaks, and will not be 
hampered by atmospheric fluctuations even in the high- 
est frequency bands. In 10 hours of integration over the 
cluster area, the sensitivity of OLIMPO in the 217 (345) 
GHz band will be « 1.5 (3.0) /iK per pixel (with a pixel 
side of 1'). These measurements will nicely complement 
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Fig. 2. Left. Simulated SZ maps of the cluster 1ES0657-556 as 
observable with the SPT at v — 150 GHz (upper panel), v — 
223 GHz (mid panel), v — 350 GHz (lower panel). A value 
M x = 20 GeV has been adopted here. Right. Simulated SZ 
maps of the cluster 1ES0657-556 as observable with the SPT 
at v — 223 GHz for three different \ masses: M x = 20 GeV 
(upper panel), 40 GeV (mid panel) and 81 GeV (lower panel). 



those of the ground based instruments, allowing an effec- 
tive separation of the different astrophysical components, 
with special sensitivity to the fluctuating background from 
unresolved galaxies. An exciting perspective is the up- 
grade of a space-borne telescope with a Fourier Transform 
Spectrometer. This can provide the combination of high 
throughput and spectral resolution needed to compare the 
SZ t h null to the surrounding frequencies, thus providing 
convincing evidence to distinguish between the different 
signal components (see FigH]). A first test can be carried 
out already with the OLIMPO experiment, preparing the 
way for a future satellite mission. 

3. Conclusions 

The SZdm effect is an inevitable consequence of the pres- 
ence and of the nature of DM in large-scale structures. Its 
analysis in the special case of 1ES0657-556 can provide a 
direct physical probe for the presence and for the nature 
of DM in cosmic structures. 

How these observations compare to other possible DM sig- 
nal detections? In general, both gamma-ray and radio ob- 
servations of DM annihilation are quite powerful to set 



constraints on the neutralino mass and composition (see, 
e.g. Colafrancesco et al. 2006, 2007). However, in the case 
of 1ES0657-556 the expected gamma-ray emission associ- 
ated to the DM clumps is too low ( ^ 1 count vs. ~ 10 
background counts at E > 1 GeV) and cannot be re- 
solved by GLAST from other possible sources of gamma- 
ray emission, both from the cluster 1ES0657-556 and from 
AGNs in the field. In addition, the GLAST spatial resolu- 
tion (~ 9 — 18 arcmin at 10 and 1 GeV, respectively) can- 
not provide any clear spatial separation between the DM 
gamma-ray signals (expected to be concentrated at the 
DM clumps) and other possible gamma-ray signals origi- 
nating within the atmosphere of 1ES0657-556 . Radio tele- 
scopes have, in principle, excellent resolution and sensitiv- 
ity to probe the different spectra and brightness distribu- 
tion of the DM-induced synchrotron emission. For the sake 
of illustration, we evaluated that the DM-induced syn- 
chrotron emission from the East DM clump is ~ 3 — 10 mJy 
(for a smooth or smooth plus 50% mass dumpiness NFW 
DM profile, M x = 40 GeV, {av) = 4.7 • lO^cra 3 ^ 1 
model used by Colafrancesco et al. 2006 for Coma, with a 
1 fiG magnetic field) at v — 100 MHz, still marginally de- 
tectable by LOFAR. Theoretical uncertainties associated 
to the amplitude of the magnetic field in the DM clumps 
of 1ES0657-556 render, however, the prediction of the ex- 
pected signals quite uncertain. 

In such a context, the possible detection of the SZjjm ef- 
fect will provide an important complementary, and maybe 
unique, probe of the nature of DM. 

Acknowledgements. The author thanks the Referee for useful 
comments and suggestions. 

References 

Barrena R. et al., 2002, Afe A, 386, 816 
Clowe, et al. 2006, preprint [astro-ph/0608407| 
Colafrancesco, S. & Mele, B. 2001, ApJ, 562, 24 
Colafrancesco, S. et al. 2003, A&A, 397, 27 
Colafrancesco, S. 2004, A&A, 422, L23 

Colafrancesco, S. Profumo, S. & Ullio, P. 2006, A&A, 455, 21 
Colafrancesco, S. Profumo, S. & Ullio, P. 2007, PRD, in press 
Giovanelli, R. et al. 1998, AJ, 116, 2632 

Glenn J. et al. 1998, Proc. SPIE 3357, 326, Advanced 

Technology MMW, Radio, and Terahertz Telescopes 
Gomez P., Romer A.K., Peterson J.B. 2003, |astro-ph/0311263] 
Kosowsky, A. 2003, New Astron. Rev., 47, 939 
Lagache, G Dole, H., Puget, J.L. 2004, ApJS, 154, 112 
Lee A. et al. |http:/7 bolo. berkeley.edu/apexsz/instrument. html 
Liang, H. et al. 2000, ApJ, 544, 686 
Markevitch, M. et al. 2002, ApJ, 567, L27 
Markevitch, M. et al. 2004, ApJ, 606, 819 

Masi S., Calvo M., Conversi L., et al., 2005, in Proc. of 
the Intern. School of Physics E. Fermi, Course CLIX, 
F.Melchiorri & Y.Rephaeli Eds., SIF, IOS press, p.359. 
Navarro, J., Frenk, C. & White, S.D.M. 1997 ApJ, 490, 493 
Petrosian, V., Madejski, G. and Luli, K. 2006, ApJ, 652, 948 
Ruhl J., et al., 2004, Proc. SPIE, 5498, p. 11 
Runyan M.C. et al., 2003, New Astron. Rev., 47, 915 



